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ABSTRACT 

In spiral galaxies, we explain their non- 
Keplerian rotation curves (RCs) by means of 
a non-luminous component embedding their 
stellar-gaseous disks. Understanding the detailed 
properties of this component (labelled Dark 
Matter, DM) is one of the most pressing issues of 
Cosmology. We investigate the recent relation- 
ship (claimed by Walker et al. 2010, hereafter 
W+10) between r, the galaxy radial coordinate, 
and Vh{r), the dark halo contribution to the 
circular velocity at r, a) in the framework of 
the Universal Rotation Curve (URC) paradigm 
and directly b) by means of the kinematics of 
a large sample of DM dominated spirals. We 
find a general agreement between the W+10 
claim, the distribution of DM emerging from the 
URC and that inferred in the (low luminosity) 
objects of our sample. We show that such a 
phenomenology, linking the spiral's luminosity, 
radii and circular velocities, implies an evident 
inconsistency with (naive) predictions in the A 
Cold Dark Matter (ACDM) scenario. 

Keywords: galaxies: spiral, kinematics and dynam- 
ics; dark matter 

I. INTRODUCTION 



Rotation curves (RCs) of spiral galaxies show no 
asymptotic Keplerian behavior and fail to match the dis- 
tribution of the luminous matter. The favored explana- 
tion is the existence of an unknown unseen component, 
the Dark Matter (DM): the "luminous" components of 
galaxies such as gas and stars are distributed inside a 
spherical "dark" halo . 

Recently, [2], hereafter W+10, following the pioneer pa- 
per of [3] (M+07), found that the dark component of 
the circular velocity (and of mass) profile shows a char- 
acteristic and similar behavior both in Spirals (and in 
Milky Way dSphs) of very different luminosities. More 
precisely, they investigated a sample of 60 spirals, whose 



A thorough review on this issue can be found in pQ and down- 
loaded at www.sissa.it/ap/dmg/dmaw_presentation.html. 



RCs span radii from 1 kpc to 75 kpc, have maximum am- 
plitudes in the range 50 km s _1 ^ V max ^ 300 km s _1 
and their baryonic masses Mf, ranging between 3 xlO 8 
M and 4 x 10 11 M Q . They derived V h {r), the halo 
contributions to the circular velocity (in short the "halo 
RCs"), by subtracting the contribution of the baryonic 
matter from the circular velocity: V£(r) = V 2 (r) — V 2 (r) 
(see W+10 for details). They found a correlation between 
V h (r) and r (see W+10): 

independently of galaxy luminosity. Moreover: M ^ r ^ = 

2+ 2 12 x 10 s (j^f ■ Noticeably, at r = R 1/2 , being R 1/2 

the radius encompassing half of the stellar mass, the re- 
lation ([I]) can be investigated also in dwarf spheroidals, 
and found to hold also in these systems (see W+10 for 
details). 

These are surprising findings: in fact, when, in objects of 
very different luminosity we derive their Vh (r) (the dark 
halo- contribution to the circular velocities) and plot them 
as a function of radius, a well defined curve emerges, in- 
stead of a more likely scatter plot. On the other hand, 
the circular velocities V(r) show a very different phe- 
nomenology. From the analysis of several thousands of 
RCs, an Universal Rotation Curve (URC) emerges, i.e. 
a specific 3-dimensional surface linking a) the circular 
velocity at a normalized radius r/Ro (i.e. measured in 
units of disk length-scale) , b) the galaxy luminosity and 
c), the normalized radius (see g], [5] (PSS), [B], [7] and 
[H] (S+07) for details). 

Summarizing, in the intermediate regions of spirals, 
W+10 found: Vh(r) = F w+10 (r), independently of 
galaxy luminosity L and stellar disk length-scale Rd- 
In the URC scenario, instead, we have: Vh(r) — 
FuR.c( r I Rd, L), i.e. the halo component of the circular 
velocity is a function of disk luminosity and disk length 
scale. Let us stress that the existence of the URC implies 
that the structural parameters of the Dark and Lumi- 
nous components are strongly related, but not necessar- 
ily, through the very constraining eq. ([!]). 

In this paper, we investigate whether these two ap- 
parently different results are compatible. We also use 
the outcome to derive important information on the DM 
distribution in spirals. We will study an independent 
(larger) sample of spirals by means of reliable method 
of mass modelling. From [9] (PS95) we select 116 spi- 
ral galaxies (hereafter Sample A) whose RCs satisfy the 
following requirements: a) no presence of a prominent 
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bulge: V(0.2R opt ) < 50 km s -1 2 , or if the latter is 
missing: V(0AR opt ) < 70 km s -1 , b) low luminosity ob- 
jects: V(r) < 120 km s~ , at any radius, that implies: 
Mi > —20.2, with Mj the magnitude in the I band, c) 
the RCs have a measure either at 0.8 R opt or at 1.0 R op t- 
Let us notice that we limit ourselves to objects of low 
luminosity, in that they are DM dominated (e.g. [11] 
(PS90)) and therefore Vh(r) can be derived with great 
accuracy. Higher luminosity spirals are disk dominated 
and therefore Vh(r) strongly depends on the details and 
on the assumptions in the mass modelling. The RCs of 
Sample A have a radial range matching that of W+10 
sample max [^Rd, lkpc] ^ r ^6R D . Following W+10, 
we adopted 1 kpc as the innermost radius in that, inside 
which, the stellar disk (sometimes with a stellar bulge) 
dominates the gravitational potential making difficult the 
estimate of Vfe(r). 

Let us compare the W+10 relationship with the URC 
analogue (PSS, S+07) and with that we derive for Sam- 
ple A. From the condition of centrifugal equilibrium we 
have: 



V 2 (r) = V l l m (r) + VZ(r), 



(2) 



where the subscripts lum and h stand respectively for 
luminous and halo components. The former is the 
quadratic sum of disk, bulge and gas contributions. Here, 
we neglect the latter two because 1) the selected low lu- 
minosity objects have, for r > 1 kpc, a negligible bulge, 
2) the gaseous disk is important only for r > R op t, 
i.e. outside the region that we are considering here (see 
[12] ). For the objects of Sample A, each RC has 2-4 
independent measurements (see Table 2 of PS95) and 
therefore the same number of halo velocity data. For 
the URC, there is analytical formula for Vh(r) (see be- 
low). Hereafter, we adopt a flat cosmology with mat- 
ter density parameter VLm = 0.27 and Hubble constant 
H = 75 km s _1 Mpc -1 . 



II. THE CONTRIBUTION OF DARK MATTER 
HALOS TO THE CIRCULAR VELOCITY OF 
SPIRALS 

We derive the velocity halo contribution for the objects 
of Sample A in the following way; we assume a Freeman 
profile to describe the stellar disk surface density, then 

Vg(a:) = 14.9 P(x) U 2 (1.0) x 2 [I K - hK^, (3) 

where x — r/R opt , I n and K n are the modified 
Bessel functions. Let us also notice that: R\/2 — 



1.67/3.2 R opt = 0.52 R opt . The mass modelling of indi- 
vidual (and coadded) RCs of objects of low luminos- 
ity (as those in Sample A) is rather simple. We have: 
,3(1) ~ 0.13 and ,9(0.8) ~ 0.23 (see PS90 and PSS). More 
precisely, there is a small luminosity dependence of /3 but 
for the objects in Sample A this is irrelevant. Then, for 
any object, we derive Vh(r) by subtracting the (small) 
disk contribution V^(r), given by eq. (|3j), from the cir- 
cular velocity V 2 (r). Finally, we derive V/j(i?i/ 2 ) by lin- 
early interpolating Vh(0AR opt ) and Vh(0-6R op t)- 
The URC leads to an (analytical) Universal form for 
Vh(r), i.e. for the halo velocity component (URCH) to 
VuRc( r ) which was built by coadding the kinematics of 
thousands of galaxies (PSS, S+07). We model the URC, 
that represents the typical RC of an object of magnitude 
Mr, (or of virial halo mass M V i r (see S+07 and refer- 
ence therein) in its dark and luminous components (e.g. 
S+07) 3 ). In detail the URC DM density profile is 
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where tq is the core radius and po is the central halo 
density. It follows that: 
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From jTT] we have: 

Po = 5 x 10- 24 (r /(8.6kpc))- 1 g cm' 3 , (6) 



and 



log te 



0.66 + 0.58 log 



M„ 



10 n M P 



(7) 



The disk length-scales and masses, Rd and Md, are re- 
lated to the halo masses through (see [T5] and S+07) 



loi(^) = 0.633 + 0.37910,^ 



0.069 log 



Ml 
10 n M, 
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and 



M D = 2.3 x 1O-M ^ 3 X 1011M ^ 



3.1 



1 + [M vir /(3 x 10 11 M Q )] i 



(8) 



(9) 



By means of equations (|4])-([9]) we can derive, the URCH 
halo velocity V£ RC (r/R D ,L). 



2 Ropt is the radius encircling 83% of the light, for a Freeman ex- 
ponential disk Ropt = 3.2Ro where Ro is the disk scale-length, 
see [TOl . 



The virial mass M„,; r corresponds to the mass enclosed in a 
sphere with density 101 p c , being p c the critical density of the 
Universe. M v i r and the virial radius are related by R v i r = 
259(M Ilir /lO 12 M ) 1/3 kpc, see [13]. 
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FIG. 1: The W+10 halo velocity-radius relationship (black line) with its 1-er scatter (black dashed lines) compared with the halo 
velocity calculated at R1/2 as function of R1/2 for 1) our Sample A whose points are calculated by interpolating Vh(0AR op t) 
and Vh{0.6R op t) for each galaxy (black points) and 2) URCH profiles (red line). 



The first step is to investigate the Vh{R\/2) ~ R1/2 rela- 
tionship, that holds also for objects of different Hubble 
types. In Fig. ([I]) we plot the URCH as a red line and 
the relation for Sample A as black points. It is clear that 
both individual and URC halo velocities do correlate with 
R1/2, well matching the W+10 (M+07) relation (shown 
as black line). 

The further step is to investigate the full profile Vh(r) 
at "intermediate radii" max [^Rp, lkpc] ^ r ^ 6i?r> 
Each of the 116 spirals contributes with 2-4 independent 
kinematical measurements. In Fig. ([2| we plot Vh(r) 
obtained from the RCs of Sample A, once that the lumi- 
nous component has been removed as explained above. 
We also plot, as two red lines, the URCH velocity pro- 
files of spirals with luminosities similar to those of the 
objects of Sample A 4 , and to the great majority of the 
60 spirals of the W+10 sample. Notice that, given the 



Schcchter-likc form of the luminosity function of spirals, 
in the W+10 sample there are only an handful of big 
galaxies. 

In the low luminosity range we plot Vh{r) for the indi- 
vidual RCs of Sample A (black points) and for the cor- 
responding URCH profiles, and compare them with the 
W+10 relationship. Let us recall that we did not in- 
vestigate the individual Vh(r) of luminous spirals since 
their dark-luminous RC decomposition is quite uncertain. 
However, we have investigated the high luminosity ob- 
jects by means of the URCH, that we show in Fig. ^ 
alongside with the W+10 relationship. Summarizing, for 
spirals in halos with M V i r > 1.9 x 10 12 Af Q , we derive 
Vh(r) by means of the URCH, for halos of lower mass, 
instead, we derive Vh(r) from 1) the URCH (red lines) 
and 2) the mass modelling of the individual RCs of Sam- 
ple A (black points). In both cases, the (derived) halo 
velocities are compared with the W+10 relationship. 



The related virial masses are M vir /M Q = 1.17 x 10 11 , 7.7x 10 11 . 
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FIG. 2: The W+10 halo velocity-radius relationship (black line) with its 1-a scatter (black dashed lines) plotted with 1) the 
halo velocities of our Sample A of 116 spirals (black points) (PSS, PS90), 2) the URCH profiles (red lines), corresponding to 
objects with mass comparable with those in the W+10 sample. The error (of 20 %) in the individual determination of Vh is 
shown as an errorbar. 



III. UNIVERSAL HALO VELOCITY PROFILE 
AND IMPLICATIONS 



We find for Sample A and also from the URC, that: 
Vh(Ri/2) oc Ri'/2> ver y similar to the W+10 relation. 
This result is not a new one, it arises as effect of the 
particular dependencies of the dark and luminous struc- 
tural parameters on luminosity and disk length scale. 
We can show this by looking at the Radial Tully Fisher 
relation ([H]). From their eq. (9) one has Lj oc 
V(Ri/ 2 ) 2 55 , where L[ is the I-band luminosity, con- 
sidering that for spirals we have: R\/ 2 oc Rd oc L°-°, 
and that from the mass modelling of individual objects 
V{R 1/2 )/V h (R 1/2 ) oc V h (R 1/2 ) - 2 (see PS90, [HUH], i.e. 
smaller galaxies have a larger fraction of DM), one ar- 
rives to Vh{R\/2) oc i?2/2> W1 tb k very near to the value 
of +0.5, claimed by W+10 (see Fig. @). 
However, with the complete data from the URC shown 
in Fig. Q we realize that Vh(r) in Spirals of different 
luminosities, has not an unique profile, neither always it 



is a power law with exponent 1/2. However, at low lu- 
minosities, the halo velocity individual contributions and 
the URCH profiles arc in good agreement with the W+10 
result (see also M+07). 

At highest luminosities (see Fig. |4)) the W+10 rela- 
tionship cannot be reproduced by the URCH profiles. 
However let us stress that in these cases Vh(r) depend 
significantly (especially in the W+10 work) on the as- 
sumptions adopted in the mass modelling. Moreover, the 
URCH profiles are obtained from hundreds of RCs, while 
V h (r) in W+10 only from few (see W+10). 
As a result, in a first approximation we can claim that the 
W+10 relationship, in the radial range 1 kpc < r < 30 
kpc, is a projection, on the (r, Vh{r)) axes of a more com- 
plex VjjRc{ r I Ropt, L) relationship in which the structural 
quantities Vurc(1, L), L and Rd are related in a specific 
way. This suggests the existence of a number of scaling 
laws between disk scalenghts, disk mass, halo structural 
parameters and galaxy luminosity, which by the way have 
been also observed through the analysis of accurate mass 
models in large samples of galaxies (e.g. S+07). 
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FIG. 3: The W+10 halo velocity-radius relationship (black 
URCH profiles (red lines). 



line) with its 1-cr scatter (black dashed lines) plotted with the 



The W+10 relation URC-supported, becomes very useful 
to test whether DM halos in spirals are consistent with 
NFW density profile. In this case let us remind that (see 

my- 



Vnf w( r ) 



V 2 9(c) 
xg{cx) 



(10) 



where x — r/R V i r is the radial coordinate, 17(c) = 



[Zn(l + c) - c/(l + c)] _1 
rameter c we take 

c(M mr ) = 9.2 



and for the concentration pa- 



-0.09 



(11) 



10 12 h- 1 M Q 

(see [U |2Q] ) . The halo mass inside a radius r is given by 

9(c) 



M NFW (x) = My 



g(cx) ' 



(12) 



In Fig. Q we plot VNFw(r, M vir ) for the radial and 
halo mass ranges identical with those investigated above 
(more specifically we consider 8 masses equally log- 
spaced in the range: 7.7xlO n M vir /M ^ 1.1 xlO 13 ). 



The resulting halo velocity curves Vnfw(t) occupies a 
very well defined portion in the (r, Vh(r)) plane (blue 
lines). However, this region is mostly off that defined by 
1) the W+10 relationship, 2) the similar relationship de- 
rived in this paper by a larger number of individual RCs 
(of low luminosity objects) and 3) a specific projection 
of the URCH (given in Fig. (|2|). A disagreement of the 
NFW halo velocity profiles with the actual kinematics of 
spirals is not new (e.g. [21] ). However, here, we make 
a step further by investigating an unprecedented large 
number of low luminosity objects. We derive their Vh(r) 
profiles from disk kinematics by adopting only weak and 
reasonable assumptions. These profiles result incompat- 
ible with the NFW profiles for any value of their masses 
and their concentrations. 



IV. DISCUSSION 

We have compared the W+10 (Vh(R\/2) -R1/2) rela- 
tionship (see also M+07), with the URCH and with the 
kinematics of 116 spirals of low luminosity objects. We 
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FIG. 4: The W+10 halo velocity-radius relationship (black line) with its 1-a scatter (black dashed lines) plotted with 1) the 
halo velocities of our Sample A of 116 spirals (black points), 2) the NFW halo velocity profiles (blue lines, with the thick curves 
corresponding to objects with mass comparable with those of the individual Sample). 



find that the W+10 "unique" velocity profile is essen- 
tially in agreement with the URC paradigm of which it 
is a sort of projection. The W+10 relation alongside its 
not negligible scatter of 0.25 dex, can be identified as the 
projection of the URC on the (r, V^(r)) plane. This im- 
plies the existence of a number of scaling laws between 
disk scalenghts, disk mass, halo structural parameters 
and galaxy luminosity, observed also analyzing accurate 
mass models in large samples of galaxies. 
We claim that the distribution of dark and luminous mat- 
ter in galaxies is such that at any chosen radius r, Vh(r) 
takes, in different galaxies of different mass, almost the 
same value. The simplest explanation is that the more 
massive is the dark halo, the lesser dense and concen- 
trated it turns out to be. The net effect is that the dark 
mass inside a chosen physical radius is weakly depen- 
dent of the total halo mass. The physical explanation of 
this still lags: we realize that the dark and the luminous 
matter are distributed in a related way, which is not ob- 
vious, in view of their very different nature. It is worth 
to point out that the URCH profiles for virial masses 



M vir > 3 x 10 12 M o , are discrepant with the W+10 re- 
lationship, also considering its (quite large) scatter (see 
Fig. (|3|). Although in the W+10 sample there is an 
obvious shortage of such big objects, it would seem that 
for the about 5% top luminous objects, the halo velocity 
profile implied by W+10 does not agree with the URC 
mass model and it indicates a larger amount of DM. This 
issue will be clarified only with a bigger sample of indi- 
vidual RCs. At low masses, low circular velocities and 
low radii (1 kpc < r < 30 kpc), we find the existence 
of an "unique" Vh(r) profile whose physical meaning has 
to be understood. This is amazing in that by adopting 
the radial units, i.e. kpc, we lose any autosimilarity of 
the velocity/mass profiles, present when the radial co- 
ordinate is measured in terms of R V i r or Rjj. We can 
use the W+10 relationship as a Cosmology Test: in fact, 
when plotted in radial physical units, the NFW velocity 
profiles lie in a very small region, almost independently 
of their halo masses and concentration. Such a region, 
however, turns out to be clearly outside that indicated by 
the derived halo velocity profiles in W+10, in the URCH 
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and in our sample of individual RCs. Especially at low 
luminosities, where the NFW velocities are 0.2-0.4 dex 
higher, the observed discrepancy is very severe in that, 
being the contribution of the luminous matter quite neg- 
ligible, Vh(r) is virtually derived in a model independent 
way, and it almost coincides with the observed V(r). Of 
course, the inconsistency between RCs and NFW halo + 
disk mass model is already well known, even if it is proved 



only in a quite limited number of test cases ([22]). Here 
we have opened the possibility to investigate such a cru- 
cial issue in direct way and by means of a large number 
of objects. It is worth recalling that also within ACDM 
scenario there are several ways to modify the original 
NFW halo profiles to be compatible with observations, 
see [2"5H2"5] and references therein. 
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